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Introduction 

This contract began in April 1990 with studies on the 
effect of wingtip shape on induced drag. That work 
considered problems of determining induced drag from 
computational aerodynamic methods and examined ef- 
fects of wake roll up on the induced drag of elliptical 
and crescent wings, [1], 

The research contract was continued in April 1991 
with scope augmented to include considerations other 
than induced drag. The idea was to develop methods 
for comparing wing tip shapes on the basis of their 
effects on total drag, structural weight, and high lift 
performance. To accomplish this, work was done to 
improve methods for multidisciplinary analysis and op- 
timization of wings. 

This report describes results of research conducted 
from April 1991 through March 1992. The general ob- 
jective was to improve an existing wing optimization 
method, [2], and apply the method to specific prob- 
lems of interest. The method, while a valuable tool for 
wing tip design studies, can be applied to more general 
problems, and has been applied to some of these other 
problems during its development. 

Specific goals that were accomplished are listed be- 
low and are explained in more detail in the report. 

1. Improvements were made to the portability and 
control flow of the existing code. The major itera- 
tion loop dealing with structural design was sped 
up and an alternate approach, using the optimizer 
to do structural sizing, was studied. 

2. Analysis methods were improved in the areas of 
structural and high lift modelling. The struc- 
tural method was revised to give total wing weight 
and verified against data for particular commer- 
cial aircraft. The high lift analysis was improved 
to provide reasonable estimates of Cx mar in the 
flaps down condition. These improvements en- 
abled making wing area a design variable, where it 
had been a fixed variable in the original method. 

3. The method was applied to the design of wings for 
a Learjet. Rough studies were done to determine 
the effects of laminar flow design on wing shape. 

4. Studies on wingtip shape were begun. 

Nomenclature 

Structural Sizing: 

a Step length parameter. 

{*,} Skin thickness. 

{t JFS } Minimum skin thickness constraint. 


{/ igrf } Difference between actual and minimum 
skin thickness. 

Wing Analysis: 
c Chord. 

ChjC XXl Chh Constants for induced flap deflection fit. 
Cl Configuration lift coefficient. 

C Lmax Maximum configuration lift coefficient. 

cj Section lift coefficient. 

cjmoj Maximum allowable section lift coefficient. 

A cj Increment to section lift due to flap deflec- 

tion at zero angle of attack. 

A ci max Increment to maximum allowable section 

lift coefficient due to flap deflection. 

ci a Lift curve slope. 

c\ 6 Flap lift curve slope. 

E Young’s modulus. 

G Shear modulus. 

h Distance between points in y-z plane. 

h Local panel unit normal vector. 

t Wing thickness. 

t/c Thickness to chord ratio. 

Uoo Freestream velocity. 

v Velocity. 

x Coordinate direction, nominally along 

freestream. 

y Coordinate direction, nominally pointing 

to aircraft right. 

z Coordinate direction, nominally pointing 

up. 

a Angle of attack. 

<5 Flap deflection. 

8 g Geometric flap deflection. 

Si Flap deflection due to induced velocities. 

T Vortex strength. 

A Wing sweep. 

cr y Yield stress. 

9 Unit normal vector in direction of velocity 

induced by shed vorticity at flap edge. 

Improvements to Code 

Portability 

The code was made more portable. It runs in inter- 
active mode on the Macintosh computer and in batch 
mode on many workstations and mainframe computers, 
including Sun, Iris, DEC, and Cray Y-MP platforms. 
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Struct viral Sizing Algorithm 

A fast algorithm was developed for sizing skin thick- 
nesses. This task is an iterative process because of 
aeroelasticity: changing the skin thicknesses to meet 
stress constraints changes the elasticity of the wing, 
which in turn affects the loads and stresses. The de- 
veloped algorithm attempts to determine, in the least 
amount of time, the skin thicknesses that meet maxi- 
mum stress and minimum gauge constraints. 

Each iteration in the algorithm begins with an ini- 
tial guess of the skin thickness, { t s }\ . Structural loads 
are calculated using this guess, and a new set of thick- 
nesses are found that meet the stress constraints under 
those loads, The difference between {t SFS } 1 

and {f,}i is the skin thickness error. 

{^Jerrll = “ {^ } 1 

Perturbations that are proportional to the error distri- 
bution are considered about the initial guess. 

{M = {Ml + a 3 err } 1 

The error distribution is approximated as varying lin- 
early with the step length, a. 

{<*.„} = + « {^} 

The change in error relative to step length, | j , 
can be estimated using finite differences. The step 
length minimizing the mean square error is used to 
determine the next guess of the skin thickness. The 
mean square error can be written as follows. 

{*».,.} • = {*»«rr}l ' {^.rr}l + 

MKrrh ■ {^} + « 2 {^} ■ 

Minimizing with respect to a yields 



dt **rr \ 

da ) 

l Ja J 




The iteration steps are summarized below. 

1. Using the current guess of the skin thickness distri- 
bution, calculate loads under all structural 

design conditions. Flag the critical conditions that 
violate the maximum stress constraint. For the 
calculated loads, determine the skin thicknesses 
that meet the stress constraint, {f, FS }i> and the 
skin thickness error, {f Jerr }i* 

2. Using {* 3)2 = M fs }i f° r s ^ n Sickness dis- 

tribution, calculate loads under the critical condi- 
tions found in step 1. Determine skin thicknesses 
to meet the stress constraint, {fj FS }2, and the skin 
thickness error, {<,„„} 2 = - {* s } 2 - 


3. Calculate the change in skin thickness error. Since 
a = 1 in step 2, {%^} = {< 5. rr >2 - 
Estimate the step length, a. 



l 

da / 

1^7^ j 

f- 

[-jr} 


Evaluate the next guess for skin thickness. 

(M = {Mi + « {MMi 

The algorithm loops from step 3 to step 1 and ter- 
minates when (f Jerr }i is sufficiently small. Neglecting 
non-critical conditions in step 2 saves time by avoiding 
load cases that do not affect the required skin thick- 
nesses. The calculation of a is protected against divi- 
sion by zero errors and is prevented from excessively 
large values. 

Optimizer Directed Structural Sizing 

Using the optimizer to handle skin thicknesses as de- 
sign variables was examined as an alternate way to do 
structural sizing. This has potential advantages over 
the structural sizing algorithm described above. Opti- 
mizer directed structural sizing eliminates the internal 
iteration on skin thicknesses, leading to a simplified 
program structure, figures 1 and 2. It could also re- 
duce overall optimization time by reducing the cost of 
a function evaluation and eliminating function noise 
due to the internal iteration loop. The potential dis- 
advantage is increased optimization time due to the 
addition of design variables for the skin thicknesses. 

To compare optimizer directed structural sizing 
against internal iterative sizing, the two methods were 
run on the same problem. The primary task was to op- 
timize wing shape for minimum drag subject to fixed 
weight, low speed lift, and low speed handling con- 
straints. The specific design variables and constraints 
are listed in tables 1 and 2. 

Design variables and constraints relating to the pri- 
mary task were handled identically in the runs of both 
methods. This eliminated differences in scaling and 
step size as a factor affecting optimization efficiency. 
For auxiliary design variables and constraints associ- 
ated with the skin thickness sizing task, scaling and 
step size were selected to give fast convergence to a 
reasonable answer. 

From table 1 and the drawing in figure 3, the two 
methods appear to have found the same answer. Table 
2 indicates that constraints are closely met. The final 
value of the objective function from the two methods is 
close but not exact. This is partially due to the linearly 
varying skin thickness approximation used in the opti- 
mizer controlled method. While the internal iteration 
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element A 


wing geometry 



wing performance 
wing weight 


Figure 1: Program Structure for Iterative Skin 
Thickness Sizing. 


wing geometry 
skin thicknesses 



wing performance 
wing weight 
skin thickness error 


Figure 2: Program Structure for Optimizer Di- 
rected Skin Thickness Sizing. 



Figure 3: Optimized Wing Planforms and El- 
ement Layout. The planform found with skin 
thicknesses calculated by the optimizer, thin 
solid, is overlaid on the planform found with the 
internal iterative solution, heavy dashed. The 
layout of elements and sub elements is indicated 
to clarify the description in tables 1 and 2. 

method determines the fully stressed skin thickness at 
each panel, the optimizer directed method does this 
for select stations on the wing. Shown in figure 4, this 
results in some extra material in the wing that is not 
required for stress or minimum gauge criteria. This dis- 
advantage causes the resulting objective function from 
the optimizer controlled approach to be slightly larger 
than the result using internal iteration. 

Details of the two optimizations are shown in ta- 
ble 3. Using the optimizer in place of an internal 
iteration requires more function evaluations, mainly 
because the increased number of design variables re- 
quires additional function evaluations for gradient cal- 
culations. 

The optimizer controlled method used fewer opti- 
mization iterations, but this difference is mainly a mat- 
ter of how long the optimizer took to decide that it 
could not improve the objective function. Shown in the 
iteration histories, figure 5, the two methods converge 
to their respective optimums in a similar manner. As 
explained earlier, the final objective function values are 
slightly different because the skin thickness distribu- 
tion from the optimizer directed solution is more coarse 
than that from the internal iteration. 

In a practical sense, both optimization runs arrived 
at their optimums with the same number of iterations. 
The iterations beyond number 40 were taken to exhaust 
the optimizer’s ability to improve the objective func- 
tion. This process consumes an unreliable number of 
iterations; variations in this number can be large. This 
makes the 25% overall time savings of the optimizer 
directed method not very meaningful. The time per 
iteration is meaningful, however, because it compares 


3 










Design 

Scale 

Finite 

Internal 

Optimizer 

Variable 


Difference 

Iteration 

Directed 



Step 

Result 

Result 

BaseSemiSpan 

1.0 

1.0e-2 

14.0 

13.9 

AreaFracA 

50. 

1.0e-4 

0.535 

0.534 

AreaFracC 

100. 

1.0e-4 

3.90e-2 

3.89e-2 

TaperA 

10. 

5.0e-4 

0.643 

0.646 

TaperC 

10. 

5.0e-4 

0.459 

0.459 

SweepA 

0.1 

.02 

37.1 

36.9 

IncRl 

0.1 

.02 

1.31 

1.43 

IncT2 

0.1 

.02 

0.02 

0.11 

IncT4 

0.1 

.02 

-3.22 

-3.28 

IncR5 

0.1 

.02 

-3.03 

-3.04 

IncT5 

0.1 

.02 

-4.06 

-4.04 

TonCRl 

10. 

1.0e-3 

0.136 

0.135 

TonCT2 

10. 

1.0e-3 

0.117 

0.117 

TonCT4 

10. 

1.0e-3 

0.125 

0.124 

TonCRS 

10. 

1.0e-3 

0.122 

0.122 

TonCT5 

10. 

1.0e-3 

0.125 

0.124 

TsRl 

.01 

1.0e-3 


12.6 

TsR2 

.01 

1.0e-3 


12.0 

TsR3 

.01 

1.0e-3 


9.8 

TsR4 

.01 

1.0e-3 


2.5 

TsR5 

.01 

1.0e-3 


1.0 

TsTl 

.01 

1.0e-3 


11.9 

TsT2 

.01 

1.0e-3 


9.8 

TsT3 

.01 

1.0e-3 


2.6 

TsT4 

.01 

1.0e-3 


0.8 

TsT5 

,01 

1.0e-3 


1.0 


Table 1: Design Variables. Design variables 

used are shown with respective scaling factors, 
finite difference step sizes, and results from the 
two optimization methods. BaseSemiSpan is 
the wing semi span. AreaBYacA and AreaFracC 
are the fractions of total wing semi area allo- 
cated to elements A and C respectively. (Ele- 
ments and sub elements are marked in figure 3.) 
Taper A and Taper C are taper ratios of the ele- 
ments. AreaFracB and TaperB are determined 
from a fixed area constraint and continuity of 
chord between elements A and B. Sweep A is the 
wing quarter chord sweep. Inc, TonC, and Ts 
are, respectively, incidence, thickness to chord 
ratio, and skin thickness values for the sub ele- 
ments. Location on the sub element is given by 
the suffixed R or T: R for root, T for tip. Final 
digit suffix indicates which sub element. Values 
not given in the table are defined by linearly 
interpolating between given values. The distri- 
butions are continuous unless multiple values 
are presented for the same location. For exam- 
ple, while IncT4 and IncT5 represent incidences 
for adjacent points on the wing, their values are 
different indicating a discontinuity. 


ConstraintScale 

Limiting 

Internal 

Optimizer 



Value 

Iteration 

Directed 



(type) 

Result 

Result 

Weight 

.001 

1405.05 (e) 

1405.05 

1405.05 

CRollHL 

100 . 

-0.02434 (u) 

-0.02434 

-0.02434 

TsErrRl 

1.0 

0.0 0 ) 


-1.5e-5 

TsErrR2 

1.0 

0.0 0 ) 


-1.4e-5 

TsErrR3 

1.0 

0.0 ( 1 ) 


-2.1e-5 

TsErrR4 

1.0 

0.0 ( 1 ) 


-5.1e-6 

TsErrRS 

1.0 

0.0 ( 1 ) 


0.0 

TsErrTl 

1.0 

0.0 ( 1 ) 


-1.4e-5 

TsErrT2 

1.0 

0.0 a ) 


-2.1e-5 

TsErrT3 

1.0 

0.0 ( 1 ) 


-5.6e-6 

TsErrT4 

1.0 

0.0 ( 1 ) 


-4.3e-14 

TsErrTS 

1.0 

0.0 ( 1 ) 


0.0 


Table 2: Constraints. Constraints are shown 
with respective scaling factors, limiting values, 
and results from the two optimization meth- 
ods. The type markings for the limiting val- 
ues, shown in parentheses, indicate whether the 
limiting value was an equality constraint (e), an 
upper bound (u), or a lower bound (1). Weight 
is an index indicating the amount of bending 
material needed in the wing. CRollHL is an es- 
timate of maximum rolling moment coefficient 
due to aileron deflection at low speed. TsErrR 
and TsErrT are the difference between actual 
skin thickness and minimum constraints at sub 
element root and tip panels respectively. The 
suffixing digit indicates the affected sub ele- 
ment. (Refer to figure 3 for sub element lo- 
cations.) 


Property 

Interna] 

Optimizer 

Iteration 

Directed 


Result 

Result 

Objective Function Value 

.03484 

.03496 

Function Evaluations 

2362 

3237 

Optimization Iterations 

63 

55 

CPU Time (Cray Y-MP seconds) 

482 

360 


Table 3: Optimization Details. 
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0 . 10 . 20 . 30 . 40 , 50 . 

Spanwlse Coordinate (ft) 


Figure 4: Skin Thickness Distributions. Skin 
thicknesses sized by the optimizer are shown 
against minimal skin thicknesses required to 
meet stress and minimum gauge requirements. 
Note that the structure is virtually fully 
stressed, with deviations due only to the linear 
approximation of skin thickness between break 
points. These deviations cause the structure to 
be slightly less efficient than the structure es- 
timated by the internal solution method, and 
causes a small difference in optimal objective 
function value between the two solutions. 



0 . 10 . 20 . 30 . 40 . 50 . 60 . 70 . 

Iteration Number 


Figure 5: Iteration Histories. Both inter- 

nal iteration and optimizer directed results are 
shown. The difference in final objective func- 
tion value is due to differences in representation 
of the skin thickness distribution. (See figure 4.) 


the cost of many simple function evaluations with the 
cost of fewer evaluations possessing internal iterations. 
In this respect, the optimizer directed method was still 
better, taking 14% less time per optimization iteration 
than the internal iteration method. 

Since the optimizer directed solution uses 57% more 
function evaluations per optimizer iteration, it can be 
inferred that each function evaluation takes 55% less 
time than a function evaluation with interned iteration. 
This suggests that the per function evaluation time sav- 
ings from removing an internal iteration loop must be 
large for the optimizer directed method to be effective, 
especially if removing the loop will greatly increase the 
number of design variables. 

This study showed that it can be advantageous to 
use optimization for controlling analyses requiring it- 
erative solution. While simplifying the coding needed 
for the analysis, optimizer directed iteration can reduce 
computational time when the cost of internal iteration 
is large and the number of extra design variables asso- 
ciated with replacing the iteration is small. 


Analysis Improvements 
Structural Design 

Improvements in structural design were aimed at ac- 
curately predicting wing weight and static aeroelastic 
deformations. The resulting model with its verification 
studies are described below. 

The wing structural model includes only the box sec- 
tion between front and rear spars. All load is assumed 
to be carried in the structural box, and aeroelastic de- 
flections are based on deformations of this box. The 
box is described as a thin walled beam with vertical 
fore and aft spars and parabolic upper and lower sur- 
faces. 

For verification studies, the general box shape given 
in figure 6 was used, scaled in the vertical axis to match 
t/c . Maximum thickness was set near 35% chord; fore 
and aft spar locations and heights varied between air- 
craft. The resulting cross section provided a reasonable 
representation of actual aircraft structural box geome- 
try. 



Figure 6: Model Structural Cross Section. 

The structural box is sized in spanwise sections. 
Each section has uniform wall thickness over its cross 
section. This thickness is set by either maximum bend- 
ing stress or smeared minimum gauge. The structural 
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analysis evaluates wing loads at several potentially crit- 
ical conditions, including maneuver and vertical gust. 
The largest bending moments experienced are used for 
structural sizing. These moments are multiplied by 1.5 
to account for the ultimate load specified in FAR part 
25. Section wall thicknesses are set so that maximum 
stresses on each cross section do not exceed an allow- 
able stress, which is a fraction of the material yield 
stress. Properties of 2024-T4 aluminum are used as 
given by [3]: E = 72.4 GPa (10.5 x 10 6 psi), G = 
27.6 GPa (4.0 x 10 6 psi), <t v = 331 MPa (48 x 10 3 psi). 

The allowable stress was set at 69% of yield stress. 
This factor was chosen to provide accurate total weight 
and wing stiffness predictions for a set of commercial 
aircraft wings. In concept, the 69% factor accounts for 
weight varying proportionally with the bending weight, 
but which is not explicitly calculated. In particular, 
the correction accounts for material to resist shear and 
fatigue. 

The material outside the structural box is assumed 
to be at smeared minimum gauge. This is described 
as a skin thickness but includes material such as ribs, 
stiffeners, joints, and fasteners. In the wing box, these 
items are smeared into the walls to simplify calculation 
of bending and torsional inertias. The minimum gauge 
used is not a true minimum gauge but indicates an 
average amount of material used for each unit of wetted 
area. The minimum gauge used, 6.35 mm (0.25 in), was 
selected from correlation with weights of test wings. 

Within the structural box, sections are sized to the 
larger of the thickness required for bending strength 
or the minimum gauge thickness. This yields the box 
weight. Total wing weight is obtained by adding the 
area dependent minimum gauge weight of the structure 
outside the box. 

While testing and calibrating the analysis, the code 
was improved to account for airfoil section pitching 
moment and inertial bending relief from weight of the 
wing, fuel, and any wing mounted engines. 

Data from three different aircraft were used in the 
study. Models of the DC9-30, DC10-10, and B727-200 
were created and analyzed. The actual wing geometries 
analyzed are plotted in figure 7. Specific data on the 
geometry and flight conditions used are listed below. 


B727-200 : 


Model Data: 


Span: 

32.87 m (107.9 ft) 

Area: 

157.90 m 2 (1700 ft 2 ) 

Root chord: 

8.58 m (28.1 ft) 

Break chord 


(chord at 35% span): 

5.14 m (16.9 ft) 

Tip chord: 

2.24 m (7.3 ft) 

Sweep: 

32° 

Forward spar location: 

0.16 chord 

Aft spar location: 

0.60 chord 

Spanwise panels: 

20 

Structural Design Condition: 

Weight: 

824141 N (185267 lb) 

Altitude: 

6250 m 

Mach number: 

0.88 

Gust velocity: 

20.87 m/s (68.5 ft/s) 

Half Wing 

Results: 

estimate 

actual ratio 

wing weight 4241 kg 

4202 kg 1.01 

9350 lb 

9264 lb 

box weight 2669 kg 


5884 lb 


fuel weight 7370 kg 


16248 lb 


DC 10-10: 


Model Data: 

Span: 

47.22 m (154.9 ft) 

Area: 

358.70 m 2 (3861 ft 2 ) 

Root chord: 

13.94 m (45.7 ft) 

Break chord 


(chord at 40% span): 

7.86 m (25.8 ft) 

Tip chord: 

2.86 m (9.4 ft) 

Sweep: 

35° 

Forward spar location: 

0.15 chord 

Aft spar location: 

0.65 chord 

Engine mass: 

5731 kg (12635 lb) 


at 32.5% span 

Spanwise panels: 

20 

Structural Design Condition: 

Weight: 

1963300 N (441350 lb) 

Altitude: 

7620 m (25000 ft) 

Mach number: 

0.88 

Gust velocity: 

20.86 m/s (68.4 ft/s) 

Half Wing 

Results: 

estimate 

actual ratio 

wing weight 11175 kg 

11109 kg 1.01 

24637 lb 

24491 lb 

box weight 7990 kg 

7628 kg 1.05 

17615 lb 

16817 lb 

fuel weight 29355 kg 

33069 kg 

64716 lb 

72904 lb 
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Figure 7: The planforms described in the following aircraft models are shown from left to right: 
B727, DC10, DC9. 


DC 9-30: 


Model Data: 


Span: 

Area: 

Root chord: 

Tip chord: 

Sweep: 

Forward spar location: 
Aft spar location: 
Spanwise panels: 


28.44 m (93.3 ft) 
92.97 m 2 (1001 ft 2 ) 
5.43 m (17.8 ft) 
1.10 m (3.6 ft) 
24.5° 
0.15 chord 
0.62 chord 
20 


Structured Design Condition 1: 

Weight: 465984 N (104753 lb) 

Altitude: 7620 m (25000 ft) 

Mach number: 0.83 

Gust velocity: 20.86 m/s (68.4 ft/s) 


Half Wing Results (Condition 1): 




estimate 

actual 

ratio 

wing weight 

2381 kg 
5249 lb 

2583 kg 
5695 lb 

0.92 

box 

weight 

1506 kg 
3320 lb 

1587 kg 
3499 lb 

0.95 

fuel 

weight 

3025 kg 
6669 lb 

4762 kg 
10498 lb 



Structural Design Condition 2: 


Weight: 

522252 N (117402 lb) 

Altitude: 

7620 m (25000 ft) 

Mach number: 

0.80 

Gust velocity: 

20.86 m/s (68.4 ft/s) 

Half Wing Results 

(Condition 2): 

estimate 

actual ratio 

wing weight 2617 kg 

2583 kg 1.01 

5769 lb 

5695 lb 

box weight 1742 kg 

1587 kg 1.10 

3840 lb 

3499 lb 

fuel weight 3025 kg 

4762 kg 

6669 lb 

10498 lb 

Deflection Calculation Condition: 

Weight: 

465984 N (104753 lb) 

Altitude: 

9145 m (30000 ft) 

Mach number: 

0.75 

Wing structure: same 

as condition 2 above. 


Total weight estimates are quite good for the DC 10 
and B727. The DC9 total weight estimate varies de- 
pending on the assumed aircraft weight. The first as- 
sumed design condition is based on the original DC9- 
30 maximum takeoff weight while the second design 
condition is based on the DC9-34 maximum takeoff 
weight. Since DC9-30 and DC9-34 wings are essen- 
tially the same, there is some question as to whether 
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Spanwlse Coordinate (ft) 

Figure 8: Incidence Distributions for DC9. Cal- 
culated incidences are compared with Douglas 
estimates to check the accuracy of aeroelastic 
calculation. 

the DC9-30 wing was designed initially with aircraft 
growth in mind. The following line from Francillon, 
[4], suggests this to be the case. 


The analysis gets close results when the DC9-34 air- 
craft weight is used. 

The method slightly overpredicts wing box weight 
for the DC10 and DC9 (with DC9-34 aircraft weight). 
This is of less concern than getting accurate total 
wing weight and deflections. Wing deflections are indi- 
cated by the elastic wing incidence distributions plot- 
ted in figure 8. For comparison, figure 8 also pro- 
vides incidence distributions estimated by Douglas, [5]. 
The deflection calculation is quite good with only a 
0.4° difference in twist at the tip. This difference oc- 
curs partly because stiffening of the wing root due to 
the fuselage is not modeled. The simplification of con- 
sidering only the wing box in calculating elastic defor- 
mations also misses some stiffness contributed by ma- 
terial in the leading edge and trailing control surfaces. 
The result is a slight underprediction in wing stiffness. 
The overprediction of box weight increases the amount 
of load carrying material to improve the prediction of 
wing stiffness. 

Maximum Lift 

An accurate prediction of maximum lift is a prereq- 
uisite to selection and optimization of wing area. In 
order to free wing area as an optimization variable, 


work was done to improve the maximum lift analysis, 
particularly for predicting lift in the flaps down condi- 
tion. 

A method for predicting section c/ max with flaps 
down was implemented in the code and a method was 
developed for making corrections to c\ max due to in- 
duced camber at the flap edge. The maximum lift anal- 
ysis was verified against flight test data from DC9 and 
DC 10 aircraft. 

Basic Approach 

A critical section approach is used for predicting wing 
Ci max . The distribution of section cj is calculated and 
is compared with the local c/ max at each span wise sec- 
tion. The wing is said to be at Ci max when the lift 
coefficient on any section reaches local c/ max . 

This method gives surprisingly good results, possibly 
for the following reason: commercial transport wings 
must be designed to provide good handling qualities 
at stall. This means that stall should never begin at 
the wing tips since an asymmetric stall could cause the 
aircraft to roll out of control. In practice, because of aft 
sweep, a conventional transport wing is prone to stall 
at the wing tips, and aircraft designers modify airfoil 
sections on the inboard sections to degrade the cj maT 
of those sections and ensure stall beginning at the wing 
root. The wing is designed to stall at, or just below, 
the point where the critical section reaches its c \ max . 

The specific procedure used compares the streamwise 
ci calculated from the 3D Weissinger model with the 2D 
c/ mar of the streamwise section. No conversion is made 
to sections perpendicular to the wing sweep axis. This 
is still a topic under study, but the assumption made 
during this work was that the viscous boundary layer 
flow follows the streamwise direction more closely than 
the sweep perpendicular direction, making a compari- 
son between streamwise c/ and ci max a better indicator 
of stall than a comparison in the sweep perpendicular 
direction. 

To simulate a wing with flaps deployed, 2D section 
data is again employed. In 2D, flap deflection increases 
the ci at zero angle of attack and increases ci max . We 
note these changes with Ac/ and A C{ max respectively. 

Ac/ = Ac/(<5) 

Aci max = Ac; mal (<5) 

Ac/ and A c/ max are functions of 6 and partially reflect 
viscous effects. In translating this to 3D, incidences on 
the Weissinger model are incremented an amount, Aa, 
that gives a lift increment equivalent to that generated 
by the flap. 





From the onset of the DC-9 programme, 
the aircraft had been planned for substantial 
growth as its MGTOW was initially limited to 
less than 80,000 lb to fit within the then cur- 
rent FAA restrictions for operations by a two- 
man crew whilst it had to use JT8D turbofans 
derated from their normal take-off thrust of 
14,000 lb (6,350 kg) to 12,000 lb (5,443 kg). 
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This A a is applied to the streamwise airfoil section and 
the Weissinger model automatically handles effects due 
to sweep. Calculation of section stall is done as with 
unflapped sections except that Acj max is added to the 
section c lmax . 

When this approach is used to calculate Ci max for a 
wing with flaps deployed, a problem occurs. The crit- 
ical section appears adjacent to the flapped portion of 
the wing, on the unflapped portion. The flap induces 
upwash on this section, increasing section lift. This 
upwash has a chordwise variation that alters the effec- 
tive camber of the sections adjacent to the flap. The 
change in effective camber should increase the section 
c/max > b ut in the standard critical section analysis, the 
ci max of the geometric section is used without correc-_ 
tion for induced camber. The load calculation captures 
the extra c\ due to flap upwash, but when compared 
with uncorrected section ci maxi stall is predicted pre- 
maturely. The simple critical section approach under- 
predicts CLmax w i n g s with flaps deployed. 


Induced Camber Estimation 

To correct for this problem, a rational method for es- 
timating the induced camber is desired. From this, 
adjustments to the ci max distribution could be made 
to yield more accurate predictions of wing Cx mflX . 

The first step toward developing a method for esti- 
mating the induced effects is to identify the primary 
source of the induced velocities as the shed vorticity 
concentrated around the flap edge. Induced veloci- 
ties are proportional to the total shed vorticity, which 
is proportional to the extra vorticity carried on the 
flapped section. Induced angles are proportional to the 
shed vorticity divided by the Recalling c\ = ^ 7 , 
the induced angles are proportional to the extra c\ car- 
ried on the flapped section. 

The shed vorticity at the flap edge is actually dis- 
tributed so the lift distribution is continuous and there 
are no infinite induced velocities. Despite this, the shed 
vorticity is sufficiently concentrated that the induced 
velocities may be approximated by a radially symmet- 
ric spatial relation, as would be the case for a true 
vortex. 

From these observations, a function for estimating 
an induced flap deflection can be proposed. 


l + c k (A) + c„(f)» + c* k (f)aJ 

Here Ac/ g is the increment in 2D section c\ due to ge- 
ometric flap deflection. It indicates the strength of the 
vorticity trailing from the flap edge. C\ $ is the 2D flap 
lift curve slope. It properly converts c\ to flap deflec- 
tion. The coefficients, c/», c xr , and c ^ are fit constants 



to be determined by comparison with test calculations 
of 6i for representative geometries, x and h are spa- 
tial factors. They are components of a position vec- 
tor between three-quarter chord points of the flap edge 
section and any section, x gives the distance in the x 
direction; h give s the dis tance perpendicular to the x 
direction: h = \Jy 2 + z 1 . The dot product, 6 • xi, gives 
Si the correct sign according to the direction of induced 
velocity, 9 and the local panel normal, n. The induced 
flap angle changes from negative to positive when mov- 
ing from the flapped to unflapped wing sections. 

A set of Lin Air models were used to obtain esti- 
mates of (5,*. The models had very large aspect ratios to 
simulate “infinite” wings. Different sweep angles were 
used to determine sweep effects. Two chordwise panels 
were used, simulating the main wing section and a 25% 
chord flap. This configuration is fairly representative 
of commercial transport wings. 

For this paneling, in 2D, velocities induced at the 
control points by the vortices are 




V 2 


1 

27 re 
1 

27 re 


[f r i — 4 r 2 ] 

flri + sr*] 


Induced velocities are equated to free stream veloci- 
ties through the control points, enforcing tangent flow 
boundary conditions. 


Vi = a Uoo 

t > 2 = (a + < 5 ) t/00 


This yields a set of equations to be solved for T ; and 

r 2 - 


fi 


U a 



2fl~c(12a -f 46) 


80/3 


r 2 _ 2ttc( fa + |j) 
an Uoo 80/3 


Before proceeding further, the problem should be 
translated to 3D. Since, for now, we are dealing only in 
inviscid flow, simple sweep theory holds. The following 
transformations apply: 



a 


6 


ccos A 

Uoo COS A 
a 

cos A 

S 

cos A 


When applied to the equation for Ti and T 2 , the net 
change is that the chord length is replaced by c cos A. 


1 

27 rc cos A 


fs 


u 



u a 
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The solution for section T\ and f 2 in 3D is 


0° Swept Wing 


r i _ 27rccos A(12a + .45) 

UZ~ 80/3 

r 2 27rccos A(|a + ~S) 

U^ = 80/3 

where a and 6 sure referenced to the streamwise section. 
Lift coefficients of the panels can be evaluated using 

2r 


From this we obtain 


2w cos A / Q o _x 

^ = -3/4-0+^) 


c h 


27TCOS A 

1/4 




The Lin Air analysis provides values of c\ x and c/ 3 
from which effective a and 8 may be backed out using 
the above equations. 


2 ci x — c\ 2 

a — — 1 - 

47 r cos A 

6 _ 3^ - C), 

47r cos A 

The effective flap deflection, <5, is composed of parts 
due to geometry, 8 g , and induced velocities, <5,-. 


6 — b g + <5j 


By subtracting 6 g from the 6 calculated from Lin Air, 
the induced 5,* can be obtained. 

Three wings were modeled in Lin Air to get esti- 
mates of effective S. Sweep on these wings, shown in 
figure 9, varied from 0°, to 20°, to 40°, The spatial fit 
coefficients, c^, c XXl and chh , were selected to give close 
estimates of <5*, giving the following results. 


Ch — 13, c xx — 151, Chh — 13 


In figures 10 through 12, estimates of 6 from the spatial 
fit are compared with calculated values from Lin Air. 
The spatial fit does a good job at matching the data 
in these examples. 

To justify the use of this quasi-empirical estimate for 
5{ it must be applied to a realistic configuration. Lin 
Air was again used to get data on effective 6 , but on a 
simulated DC9 planform. This model introduced span- 
wise variation in chord and 36% chord flaps as opposed 
to the 25% chord flaps used in the previous models. 

Variation in chord is expected to affect Si somewhat. 
However the induced velocities that cause Si diminish 
very quickly on either side of the flap edge. Over the 



Figure 9: Model Planforms for Determining In- 
duced Flap Deflection. 



Figure 10: Effective Flap Deflections on 0° 

Swept Wing. Incidences calculated with Lin Air 
are compared with estimates generated using 
the proposed fit to <5; . 



Figure 11: Effective Flap Deflections on 20° 
Swept Wing. Incidences calculated with Lin Air 
are compared with estimates generated using 
the proposed fit to 5,. 
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Figure 12: Effective Flap Deflections on 40° 

Swept Wing. Incidences calculated with Lin Air 
are compared with estimates generated using 
the proposed fit to <5,. 



Figure 13: Effective Flap Deflections on Sim- 
ulated DC9 Wing. Incidences calculated with 
Lin Air are compared with estimates generated 
using the proposed fit to £ t . 

local region where the induced velocities are important, 
the variation in chord is usually small, enabling effects 
of taper to be neglected. 

The change in flap chord affects the calculation of <5 
from Lin Air. The paneling in this case is the same as 
before, except that the main element uses the first 64% 
chord and the flap element uses the last 36% chord. 
This changes the equations for v\ and V 2 , and the sub- 
sequent steps can be repeated to get the following equa- 
tion for 5. 

s = U4822^ 

4x cos A 

c\ 6 also differs from the result for 25% chord flap. 
c h = 1 .2558 1395 l7r 

When the above modification to the 2D section anal- 
ysis are used with the spatial fit to <5» developed earlier, 
a good estimate of effective 6 can be achieved as shown 
in figure 13. 



.00 .05 .10 .15 .20 .25 .30 

t/c 


Figure 14: Model for maximum section lift co- 
efficient as a function of t/c . The curves are 
estimates to c/ max for families of airfoils based 
on the DC9 and DC10 sections. 



Figure 15: Model for the increment in section 
lift due to flap deflection, Ac/, and the incre- 
ment in maximum section lift, Ac/ max . The dif- 
ference between the curves is due to different 
flap chord to wing chord ratios: 36% on DC9, 
32% on DC10. 

Comparison with Flight Test 

The final test of the method for estimating induced flap 
incidence was to use the formula for <5,- within the wing 
optimization program and attempt to match flight test 
data for the DC9 and DC10. Empirical data were used 
for c fmox > and Ac imai ■ Basic section Ci mai was es- 
timated from data provided in [6] and [7]. Increments 
to section lift, Ac/ and Ac/ max , were estimated using 
methods suggested by [8]. The models are described in 
figures 14 and 15. The resulting aircraft CL maj . is plot- 
ted as a function of flap deflection in figure 16, where 
flight test data from [9] is also provided for comparison. 

Application of Method 

With the modifications to the analysis methods listed 
above, the optimization program was used for a study 
exploring the effect of design for laminar flow on a busi- 
ness jet. This study was done to assist work in a coop- 
erative agreement between NASA and Learjet. 
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Figure 16: Maximum lift calculated by the wing 
optimization program using the <5* correction 
are compared to flight test data. The variation 
with flap deflection is shown for DC 9 and DC 10 
aircraft. Configurations examined are with slats 
retracted. 


In the study, wings were optimized for minimum 
cruise drag subject to fixed weight and minimum safe 
speed constraints. The specific design conditions are 
listed below. 


Cruise Design Condition: 

Weight: 73371 N (16494 lb) 

Altitude: 10668 m (35000 ft) 

Mach number: 0.82 

Structural Design Condition: 


Weight: 

Altitude: 

Mach number: 
Maneuver load factor: 
Gust velocity: 


73371 N (16494 lb) 
7620 m (25000 ft) 
0.88 
2.5 

20.0 m/s (65.6 ft/s) 


Minimum Speed Requirements: 

Clean: 73.7 m/s (143 kt) 

Flaps down: 64.5 m/s (125 kt) 


The aircraft is expected to be designed for a cruise 
Mach number of 0.85; however the compressibility drag 
routines implemented in the code are for old peaky air- 
foil sections and it is expected that supercritical airfoil 
sections will reduce divergence Mach number by 0.03. 
This explains the choice of 0.82 cruise Mach number 
for the study. 


To simulate effects of laminar flow, the value of skin 
friction coefficient was reduced in the parasite drag 
model. Skin friction coefficient was reduced 20% as an 
estimate for a 40% run of laminar flow. In exchange 
for skin friction reduction, wings designed for laminar 
flow were given sweep constraints. These constraints 
represent difficulties in achieving laminar flow at high 
sweep angles because of cross flow instability. Lami- 
nar flow is assumed only on the outer wing elements. 
Tbrbulent flow, and its higher c/, are assumed on the 
inboard wing elements. Laminar flow is not expected 
on the inboard elements because of contamination from 
the fuselage boundary layer. To minimize any loss from 


this, no constraint is placed on inboard section sweep. 
To obtain laminar flow on the outboard wing panels, it 
is expected that a boundary layer fence will be placed 
between the turbulent inboard and laminar outboard 
sections. 

The results of the study are shown in figure 17 and 
table 4. This study was the first to use the method 
with wing area as a design variable. The previously 
mentioned improvements to structural and maximum 
lift analyses were crucial to making proper optimiza- 
tions of wing area, and the results of this study suggest 
the necessity of such optimization. 

The laminar flow designs have more wing area than 
their conventional counterparts. Designs D, E, and F 
have 8%, 16%, and 17% more area than designs A, B, 
and C, respectively. Compared with the turbulent flow 
designs of similar weight, the laminar flow wings have 
less parasite drag, and more compressibility drag be- 
cause of the sweep constraint. Since the cost of adding 
wing area is small for laminar flow wings, the optimizer 
increases wing area to reduce c/ and hence compress- 
ibility drag. Had this been a fixed area study and the 
laminar flow wings were forced to have the same area 
as the conventional ones, the laminar flow design strat- 
egy would have been unfairly disadvantaged, stuck with 
less than optimal wing area. 

One fixed parameter that may have an important ef- 
fect on this study is cruise altitude. A fixed altitude 
constraint is reasonable to assume as it is either set by 
operational considerations or by engine performance. 
Comparing the drag of wings at different altitudes is 
of limited utility since engine performance is very sen- 
sitive to altitude. The 35000 ft cruise altitude used 
in this study may be on the low end for desired busi- 
ness jet operations, but it probably serves as a good 
altitude to compare laminar and turbulent flow de- 
signs. As wing weight increases, designs presumably 
move toward aerodynamic optimums: the area on tur- 
bulent flow wings decreases, indicating minimum drag 
at higher Cl ; the area on laminar flow wings increases, 
indicating minimum drag at lower Cl • The perfor- 
mance of laminar flow wings relative to conventional 
wings will likely improve at lower altitude and dimin- 
ish at higher altitude. 

The study was done initially with a 20° constraint on 
wing sweep for laminar flow. This was perceived as a 
sweep angle which could be reasonably expected to give 
laminar flow, but as indicated in figure 17, laminar flow 
design provides no advantage over conventional design 
with this sweep constraint. Indicated in table 4, the 
sweep constraint creates a large compressibility drag 
penalty, more than three times the compressibility drag 
on turbulent flow wings of similar weight. Wing area 
is increased to compensate, but this reduces structural 
efficiency, leading to shorter spans and higher induced 
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Turbulent Flow Cases Laminar Flow Cases 



Design A: 

Design B: 

Design C: 

Design D: 

Design E: 

Design F: 

Design G: 

Design H: 

Design I: 

BaseSemiArea 

9.86 

9.79 

9.75 

10.65 

11.31 

11.42 

9.96 

9.97 

9.98 

BaseSemiSpan 

5.68 

6.41 

6.98 

5.24 

5.72 

6.26 

5.62 

6.36 

7.00 

AreaFracA 

0.388 

0.388 

0.388 

0.387 

0.388 

0.389 

0.388 

0.389 

0.389 

AreaFracB 

0.331 

0.332 

0.332 

0.331 

0.332 

0.332 

0.329 

0.331 

0.331 

TaperA 

0.837 

0.841 

0.833 

0.870 

0.888 

0.892 

0.838 

0.920 

0.909 

Sweep A 

37.1 

37.4 

37.6 

36.1 

36.8 

36.9 

36.7 

37.2 

37.4 

IncRl 

0.63 

0.46 

0.28 

1.74 

1.32 

1.33 

1.10 

1.14 

0.98 

IncT2 

0.51 

0.10 

-0.33 

-2.70 

-1.63 

-1.76 

-1.24 

-0.97 

-1.34 

IncT3 

-1.88 

-1.45 

-0.92 

-4.21 

-3.00 

-2.93 

-2.47 

-2.73 

-2.23 

IncT4 

-4.25 

-3.66 

-3.13 

-4.59 

-3.32 

-3.12 

-4.37 

-3.33 

-2.85 

TonCRl 

0.138 

0.135 

0.133 

0.140 

0.137 

0.136 

0.138 

0.130 

0.129 

TonCT2 

0.112 

0.114 

0.115 

0.096 

0.088 

0.087 

0.109 

0.106 

0.106 

TonCT3 

0.117 

0.117 

0.115 

0.095 

0.090 

0.089 

0.109 

0.110 

0.109 

TonCT4 

0.127 

0.125 

0.123 

0.112 

0.103 

0.102 

0.135 

0.127 

0.125 

ftot 

0.332 

0.306 

0.290 

0.374 

0.337 

0.318 

0.313 

0.287 

0.270 

finv 

0.114 

0.091 

0.076 

0.136 

0.110 

0.093 

0.114 

0.092 

0.076 

fvis 

0.207 

0.206 

0.205 

0.186 

0.194 

0.195 

0.177 

0.176 

0.176 

fcpr 

0.010 

0.010 

0.009 

0.053 

0.033 

0.030 

0.022 

0.018 

0.017 


Table 4: Comparison of Laminar and Turbulent Flow Wing Designs. Optimized results for all 
design variables are shown. Tbe geometry description is similar to that used for the study in 
the Optimizer Directed Structural Sizing subsection. The element layout is the same described in 
figure 3, except only four sub elements are used: element A contains both sub elements 1 and 2, 
elements B and C contain sub elements 3 and 4 respectively. BaseSemiSpan is the wing semi span in 
meters. BaseSemiArea is the wing semi area in square meters. AreaFracA and AreaFracB are the 
fractions of total wing semi area allocated to elements A and B. AreaFracC is implicitly allocated any 
remaining wing area. TaperA is the taper ratios of element A. TaperB and TaperC are determined 
fro m continuity of chord between elements. Sweep A is the wing quarter chord sweep for turbulent 
flow designs (A to C); it is the inboard section sweep for laminar flow designs (D to I). Inc and 
TonC are incidence and thickness to chord ratio values for the sub elements. Location on the sub 
element is given by the suffixed R or T: E for root, T for tip. Final digit suffix indicates which sub 
element. Values not given in the table are defined by linearly interpolating between given values. 
Equivalent drag areas are given so designs may be compared on particular drag components. The 
posted values, ftot, finv, fvis, and fcpr, correspond to total, induced, parasite, and compressibility 
drag areas in square meters. 


drag. 

The sensitivity of drag with wing sweep was esti- 
mated for laminar flow wings by optimizing with a 30° 
sweep constraint. A reduction in drag relative to con- 
ventional designs is indicated in figure 17. It is not 
clear whether substantial amounts of laminar flow can 
be achieved on a 30° swept wing. Interpolating be- 
tween designs C, F, and I, it appears that laminar flow 
must be demonstrated on wings with at least 26° sweep 
before laminar flow design becomes useful for this par- 
ticular mission. 
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Figure 17: Comparison of Laminar and Turbulent Flow Wing Designs. 
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Figure 18: Optimized Planform with Wingtip 
Sweep as a Design Variable. 
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Figure 19: Section Lift Distributions on Opti- 
mized Wingtip Sweep Planform. 

Current Status 

As a first step toward making comparisons of different 
wing tip devices, an optimization was run with sweep 
of the tip element as a design variable. The resulting 
wing, and its distributions of lift and c/ are plotted 
in figures 18 through 20. The wing tip is planar with 
a sweep of 54° . The high sweep means the tips have 
no compressibility drag and can have high lift coeffi- 
cients in cruise. Since lift curve slope varies like cos A, 
the wing tips load more slowly with angle of attack 
than the main wing. This gives section cj’s compara- 
ble to the main wing at low speed, and leaves the tips 
relatively unloaded in the structural sizing load case. 
The optimizer found several advantages to sweeping 
the wing tips. 

While the direction for improvement is good, the ex- 
tent taken in this example is not correct because of 
effects not modeled in the analysis methods. The most 
obvious is that section c/ maa . has not been adjusted for 
sweep. 

This is not simply an oversight; there does not seem 
to be a definitive method for handling sweep effects 
on maximum lift. While some references take ci max as 
varying like cos A, [8], experiments suggest that higher 
values of c/ mar , approaching those of unswept wings, 
can be achieved, [10, 11]. 
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Figure 20: Section Lift Coefficient Distributions 
on Optimized Wingtip Sweep Planform. 

Additional work on CL max prediction is expected as 
this work continues. The current goal is to find or de- 
velop a method that accurately predicts the effects of 
sweep on section c\ max . Additional work on the struc- 
tural model is also expected, primarily on implement- 
ing an estimate for shear material, which may be im- 
portant in the analysis of swept wingtips and business 
jet wings. When these further analysis modifications 
are finished, it is hoped that practical studies of wingtip 
devices and optimal planform shapes can be made. 
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